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initio  calculation  of  the  ground  state  of  the  r^t- ’molecule  is 
performed  and  is  the  major  contribution  of  this  work.  The  implica- 
tions of  this  calculation  on  the  breakdown  problem  and  on  other 
problems  are  discussed.  /\(This  calculation  and  its  Implications 
form  the  doctoral  Jthesis!  of  Jack  K.  Lloyd,  Jr.) 


DIELECTRIC  BREAKDOWN  IN  SOLIDS 
Paul  P.  Budenstein 
Physics  Department 

Auburn  University,  Auburn,  Alabama  36830 

The  goal  of  this  study  on  dielectric  breakdown  in  solids  was  to  develop 
a theory  of  breakdown  that  was  consistent  with  experimental  evidence.  The 
initial  idea  for  the  theory  came  during  the  summer  of  1973  when  the  princi- 
pal investigator  was  reviewing  breakdown  phenomena  under  nanosecond  excita- 
tion(l).  The  proposed  theory  assumes  that  breakdown  is  initiated  when  a 
localized  charge  density  reaches  a critical  value.  Then  chemical  bonds  are 
disrupted  and  atoms  of  the  solid  are  fTeed  to  form  a high-temperature,  high- 
pressure  gas.  Electrons  are  assumed  to  be  released  during  the  bond  breakage 
and  later  in  reactions  of  the  gaseous  atoms.  These  cause  further  breakdown 
of  the  chemical  bonds.  Ultimately,  through  continuation  of  the  above  reac- 
tions and,  possibly,  through  crack  propagation,  a gaseous  channel  grows  in 
the  dielectric  until  it  bridges  the  electrodes.  Only  then  does  the  high 
current  of  breakdown  occur.  Joule  heating  subsequent  to  the  formation  of  the 
channel  can  cause  further  erosion  of  the  solid. 

If  the  above  is  essentially  correct,  then  the  fundamental  issue  of 
dielectric  breakdown  of  solids  is  to  properly  characterize  the  conversion 
of  the  solid  to  the  gas.  During  the  academic  year  of  1973-1974,  neither 
the  principal  investigator  nor  his  graduate  students  were  able  to  work  on 
the  breakdown  problem.  Thus  an  extension  of  time  was  requested  to  allow  the 
work  to  be  carried  on  during  the  1974-1975  year.  This  request  was  generously 
granted.  In  the  summer  of  1974,  the  principal  investigator  again  worked  on 
the  breakdown  problem  at  the  Guidance  and  Control  Laboratory  of  the  U.  S. 

Army  Missile  Command,  Redstone  Arsenal,  this  time  designing  breakdown  experi- 
ments employing  a megavolt,  nanosecond  pulser  with  thick  solid  samples  and 
also  pursuing  the  implications  of  his  proposed  breakdown  model.  The  latter 
have  been  summarized  in  a technical  report  (2)  and  in  two  conference  presen- 
tations (3,4).  In  the  fall  of  1974,  two  graduate  students  (Jack  Lloyd  and 
James  Knaur)  undertook  portions  of  the  breakdown  study. 

Lloyd,  a doctoral  candidate  with  considerable  experience  with  computers 
and  a fairly  good  background  in  solid  state  theory,  was  given  support  funds 
under  this  contract  to  consider  the  theoretical  question  concerning  the 
magnitude  of  the  excess  charge  density  required  for  bond  disruption  and 


2 


the  energies ( of  the  released  atoms.  Knaur,  using  the  megavolt, 
nanosecond  pulser  o£  the  Guidance  and  Control  Laboratories,  set  out  to 
make  more  detailed  observations  of  the  electrical  waveforms  and  light 
emissions  than  was  previously  possible.  This  work  was  done  partially 
to  provide  certain  engineering  data  needed  by  the  Missile  Command  and 
partially  to  provide  information  pertinent  to  the  breakdown  mechanism. 

Basic  support  for  Knaur' s work  has  come  from  the  Guidance  and  Control 
Laboratory,  but  a small  portion  of  the  funds  of  the  current  contract 
have  been  used  to  purchase  specimen  materials  and,  since  Knaur' s work 
is  the  basis  for  his  M.  S.  thesis,  Auburn  University  has  contributed 
machine  shop  time  to  fabricate  samples  and  fixtures.  Also  two  spectro- 
meters and  associated  equipment  have  been  loaned  by  Auburn  University 
for  use  at  the  test  facility  at  Redstone  Arsenal.  During  the  summer 
of  1975  the  principal  investigator  was  provided  with  funds  to  work  20 
days  with  Knaur  at  Redstone  Arsenal. 

A third  graduate  student,  Aradhana  Baruah,  was  added  to  the  break- 
down program  in  July  of  1975  and  given  some  financial  aid  from  whis  con- 
tract. Her  problem,  which  is  to  be  her  doctoral  dissertation,  is  to  make 
detailed  calculations  of  the  growth  of  the  gaseous  region  using  the  chemical 
kinetics  model  described  in  Ref.  2.  The  chemical  kinetics  model  involves 
generation  of  the  gas  from  the  solid,  reactions  in  the  high  temperature 
gas,  and  heat  conduction  from  the  gas  to  the  solid.  From  the  model  one 
obtains  such  information  as  the  rate  of  growth  of  the  gaseous  region,  its 
changing  composition,  and  its  temperature  and  pressure  vs  time.  From  the 
temperature,  the  degree  of  excitation  of  the  gas  can  be  inferred  and  hence 
its  light  emission  spectrum.  One  of  the  major  goals  of  these  calculations 
is  to  delineate  the  material  parameters  that  are  important  in  the  breakdown 
process.  It  is  hoped  that  the  calculations  will  ultimately  lead  to  expres- 
sions for  the  thickness  dependence  and  temperature  dependence  of  the  break- 
down voltage. 

The  status  of  the  three  efforts  will  be  summarized  here.  Lloyd's 
work  constitutes  the  major  portion  of  this  contract  and  is  nearing  comple- 
tion. It  will  be  his  doctoral  dissertation  and  should  provide  the  basis 
for  a substantial  journal  article.  His  thesis  and  the  article  will  be 
submitted  as  the  major  documentation  for  this  contract.  Calculations  for 
this  work  are  just  about  completed  and  work  on  the  final  document  is  in 
progress. 
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Lloyd  first  attempted  to  modify  the  traditional  approaches  of  band 
theory  to  analyze  the  effect  of  localized  charge  distributions  on  the 
realignment  of  energy  levels.  However,  the  problem  soon  became  too  com- 
plex and  an  alternate  approach  was  sought.  We  decided  the  basic  question 
was  the  alteration  of  chemical  bonds  by  the  addition  of  extra  electrons 
and  that  this  question  could  be  attacked  by  examining  bonding  in  simple 

molecular  systems.  The  simplest  molecule  is  the  1L  molecule.  Theoretical 

• • + “ 
solutions  previously  existed  for  the  H2  ion,  the  II2  neutral  molecule, 

and  the  H 2 molecular  ion;  all  of  these  are  stable  in  the  ground  state. 

Lloyd  has  performed  an  ab  initio  calculation  of  the  ground  state  of  the 
H2  molecular  ion.  A valence-bond  method  similar  to  that  used  by  Fischer- 
Hj aimers  for  the  calculation  of  the  ground  state  of  the  H2~  ion  was  employed. 
Only  Is  electronic  basis  functions  were  used,  but  with  a choice  of  differ- 
ent exponents  in  the  orbital  wavefunctions  for  greater  accuracy.  Briefly, 
the  method  entails  the  use  of  the  basis  functions  to  form  a wave  function 
having  the  proper  symmetry  and  then  computing  the  state  energy  from  the 
expression  F /ip*  HUUV 

[ ij)*  LjJd\/ 

The  internuclear  distance  appears  as  a parameter  in  the  wave  function  and 
also  in  the  energy  expression. 

Lloyd  is  making  his  final  calculations  now  and  is  plotting  his  results. 
With  several  terms  still  to  be  included,  he  obtains  a non-bonding  state  at 
the  internuclear  separation  corresponding  to  the  equilibrium  position  of 
the  neutral  H2  molecule.  The  energy  should  asymptotically  approach  that 
for  two  single  H ions  in  the  limit  of  infinite  separation.  The  height  of 
the  curve  at  the  point  of  maximum  stability  for  the  neutral  system  will 
give  some  indication  of  the  degree  of  instability  of  the  neutral  system 
with  respect  to  the  addition  of  the  two  extra  electrons.  The  curve  will 
also  allow  estimation  of  the  energy  liberated  per  particle  upon  bond 
disruption. 

Kinetics  of  cavity  growth  were  formulted  in  Ref.  2 (see  Appendix  A, 
taken  from  Ref.  2),  but  the  details  of  the  model  were  not  quantified.  The 
first  task  toward  quantification  has  been  performed  by  Baruah.  She  has 
made  a literature  study  of  very  fast  reactions  that  take  place  at  high 
temperatures  and  at  **'’  ,h  pressures  to  determine  reasonable  parameters  for 
the  rate  equations,  an  appropriate  equation  of  state  for  the  cavity  gas, 
and  an  internal  energy  function  for  the  cavity  gas.  Equations  4a  to  4d 
of  Appendix  A describe  the  rates  of  change  of  each  of  the  constituents 
of  the  gas  in  the  highly  idealized  model.  Only  four  constituents  are 
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assumed  to  exist:  neutral  atoms,  positive  ions,  negative  ions,  and  free 

electrons.  Once  the  atoms  and  ions  have  been  released  from  the  wall,  the 
reactions  are  taken  to  be  ionization  of  neutrals,  recombination,  charge 
transfer  from  a negative  ion  to  a positive  ion,  charge  exchange  between 
two  neutrals,  electron  attachment  to  a neutral,  and  electron  detachment 
from  a negative  ion.  Three  types  of  kinetic  coefficients  occur:  wall 

release  coefficients  (G{+,  q(  , 0^Q,  0(  ) that  describe  the  portion  of 
wall  particles  released  as  positive  ions,  negative  ions,  neutrals  and 
electrons;  first  order  kinetic  coefficients  ( (3  ^ and  ) which  des- 
cribe rates  depending  ‘only  on  a single  constituent's  concentration;  and 
second  order  coefficients  ( V-?*  X3’  an<*  Y4)  which  describe  rates 


depending  upon  two  concentrations. 

We  assume  that  the  wall  coefficients  are  temperature -independent  and 
are  determined  by  the  reaction  details  when  bond  disruption  occurs.  The 
four  wall  coefficients  are  not  independent  since  they  obey  the  constraints 
of  charge  neutrality  ($  + =Q(  +C()  and  particle  conservation 

Cl  = 0/+  + C(_  + 0(o). 

The  0's  and  Y's,  however,  are  assumed  to  depend  on  the  gas 

temperature  in  the  manner  typical  of  chemical  reactions,  regardless  of 

the  kinetic  order  of  the  reactions:  (•*.;  /,  -t 

a»  L//<M 

(3.  = A;  e 


£iAr  y = A' 


1 1 

where  A^  and  A^  are  descriptive  of  collision  frequencies  and  and  E^,  are 

activation  energies.  In  the  breakdown  kinetics  problem,  very  high  rates 

of  reaction  are  required  because  the  time  scale  of  growth  of  a conducting 

channel  is  10  ^ to  10  ^ sec.  The  required  reaction  rate  seems  to  bo 

somewhat  faster  than  that  of  explosive  reactions.  For  explosive  reactions, 

such  as  the  detonation  of  TNT  or  EDNA,  the  kinetics  are  well-described 

by  Eyring's  surface  corrosion  theory  in  terms  of  a kinetic  coefficient  k 

r 


kr  = A'  T2  exp  C-AH*/kT) 

whore  A',  T^  and  MI*  are  constants  for  a given  explosive  system.  This 
equation  is  of  the  Ahrennius  typo*  Tho  kinetic  equation  defining  k^  is 

dn/dt  = k (4fTr^/s),  whore  41Tr“/s  is  the 
number  of  molecules  on  the  surface  of  a reacting  grain  of  radius  r and  dn/dt 
describes  the  change  in  this  number  per  unit  time.  Thus  k can  be  thought 
of  as  the  fractional  change  in  the  number  of  surface  molecules  per  unit  time. 
For  explosive  reactions,  k is  about  10*Vsec. 
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The  equation  of  state  for  a highly  dense,  high  temperature  gas  is 
a matter  of  conjecture.  In  explosive  reactions,  best  results  seem  to  be 
obtained  by  using  equations  that  are  perturbations  of  the  ideal  gas  law. 
The  gas  that  we  wish  to  describe  is  both  more  dense  and  at  a higher 
temperature  than  that  of  high  power  chemical  explosions.  Hence,  it  is 
plausible  to  use  the  ideal  gas  law  as  the  equation  of  state  as  a start. 
Later,  when  a more  sophisticated  kinetics  model  is  developed,  the  choice 
of  equation  of  state  will  have  to  be  re-examined. 

An  internal  energy  function  is  required  because  the  kinetic  coef- 
ficients are  temperature  dependent  and  the  temperature  depends  both  on 
the  rate  at  which  energy  is  added  during  dissociation  and  the  heat  trans- 
fer to  the  surroundings.  In  addition  to  translational  kinetic  energy  for 
each  specie,  it  is  necessary  to  consider  electronic  excitation  of  the 
neutral  atoms,  positive  ions  and  negative  ions.  The  internal  energy  may 
be  written 


Z £ic~ei/k'\ 
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(n+  + n_  + nQ  + n ) (4 T?r 73)  (1  • 5 kT)  + 
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The  term  with  the  factor  (1.5  kT)  describes  translational  kinetic  energy  of 
neutrals,  negative  and  positive  ions,  and  electrons  (assumed  to  be  in  thermal 
equilibrium)  and  the  three  summations  describe  the  contributions  of  electronic 
excitations  to  the  internal  energies  of  neutrals,  positive  ions  and  negative 
ions.  The  remaining  terms  are  descriptive  of  the  heats  of  reaction. 

This  report  is  in  the  nature  of  an  interim  report.  The  main  report  for 
the  contract  will  be  Lloyd's  thesis  and  article  for  publication.  It  is 
hoped  that  these  will  bo  ready  by  tho  end  of  this  year.  We  appreciate  the 
support  of  the  Army  Research  Office  and  hope  that  the  current  effort  makes 
a meaningful  contribution  to  the  understanding  of  dielectric  breakdown  in 
solids.  We  are  particularly  grateful  for  the  continuing  encouragement 
received  from  members  of  the  Guidance  and  Control  Laboratory,  U.  S.  Army 
Missile  Command,  Redstone  Arsenal. 
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2.  Proposed  Model  of  Dielectric  Breakdown  in  Solids 
a.  Concepts  of  the  Modal 


(1)  Critical  Charge  Density.  The  basic  assumption 
made  is  that  breakdown  occurs  when  the  density  of  excess  charge  in  an 
insulator  in  a local  region  containing  many  atoms  attains  a critical 
value.  The  source  of  this  excess  charge  may  be  field  emission  from  an 
electrode,  charge  multiplication  by  an  avalanche  process,  charge  injec- 
tion at  the  surface  of  an  ionized  void,  electron  injection  from  an 
external  accelerator,  electrons  produced  by  ionization  processes  associa- 
ted with  absorption  of  laser  light,  or  some  other  means. 


( 2)  Disruption  of  Chemical  Bonds  by  Excess  Charges.  The 
effect  of  a charge  density  greater  than  critical  on  the  solid  is  assumed 
to  be  disruption  of  chemical  bonds.  (Bowden  and  McLaren  [25],  in 

3 

attempting  to  explain  how  a relatively  small  electric  field  (4.5  X 10 
v/m)  could  initiate  an  explosive  reaction  in  a single  crystal  of  silver 
azide,  hypothesized  that  a high  local  electron  density  led  to  the  decom- 
position of  the  crystal.  Subsequent  work  on  silver  azide  has  neither 
confirmed  or  disproven  this  hypothesis  [26].) 


The  mechanism  of  bond  disruption  by  an  excess  electron  density  can 
be  visualized  by  considering  the  hydrogen  molecule  and_its  ions.  An 
electron  can  be  attached  to  U9  to  form  a stable  ion,  H,,.  The  three- 

electron,  two-nucleus  system  (Hj)  will  have  a different  set  of  electronic 
energy  levels  than  the  two-electron,  two-nucleus  system  (H2) [ 27 , 28] . 

(The  bond  length  of  ^ is  0.746  A,  while  it  is  estimated  to  be  1.65  A 
in  K,.  The  H,  molecule  is  probably  best  thought  of  with  the  two  elec- 
trons shared  by  the  two  nuclei  on  an  equal  basis,  while  the  molecule 
is  more  like  an  H weekly  coupled  to  an  H .)  The  addition  of  still 
another  electron  to  th-4.  system  [ + e -»  H,  (unstable)]  leads  to  a non- 
bonding situation.  Thio  nonbonding  may  occur  with  a variety  of  end 
products  as  illustrated  in  the  following  chemical  reactions: 


Il"“  ->  H + H + 2e 


H2_  ->  H2  + 2e 


H~"  -*  H + H"  + e 
H”  -» if  + if 
11,"  H+  + H+  + 4e 


H“~  ->  H“  + e 
H“"  ->  H*  + 3e 


6 


y 


l 


H“"  -*  H+  + H"  + 2e 
H“~  -»  H+  + H + 3e. 

The  reactions  that  dominate  should  be  the  ones  most  favorable  energeti- 
cally. These  cannot  be  determined  readily  from  calculations  because 
of  the  number  of  particles  involved.  Note  that  it  is  possible  to  end 
up  with  an  excess  of  electrons;  i.e.,  more  than  1 electron  may  be  pro- 
duced in  the  dissociation  of  , whereas  only  one  electron  was  added 
to  to  form  the  unstable  . 

(3)  Chain  Reaction.  The  third  ingredient  believed 
essential  to  the  breakdown  process  is  the  chain  reaction.  A well-known 
example  of  a chain  reaction  is  the  explosive  co-  bination  of  and  Cl^ 

gases  [29].  In  the  absence  of  light  or  an  electric  spark,  this  mixture 
is  stable  at  STP.  However,  a flash  of  ultraviolet  light  or  a spark  will 
cause  the  formation  of  HC1  with  great  rapidity,  the  reaction  being 
strongly  exothermic.  The  reaction  scheme  may  be  represented  as  follows: 

H„  + light  H + H 

4. 

H + Cl2  ->  HC1  + Cl 
Cl  + H2  ->  HC1  + H 
H + Cl  ->  HC1  + light. 

The  reaction  products  in  each  case  have  the  potential  for  causing  further 
reactions;  the  total  energy  released  is  many  times  the  initiating  energy. 
The  rate  of  reaction  depends  upon  the  concentrations  of  the  reacting  com- 
ponents, temperature,  and  pressure.  These,  in  turn,  depend  on  the  pro- 
cesses of  diffusion,  loss  through  boundaries,  heat  dissipation,  wall 
effects,  competing  reactions,  etc. 

In  the  breakdown  problem,  we  assume  that  an  electron  density  is 
created  in  a local  region  sufficient  to  disrupt  chemical  bonds.  Atoms 
are  thus  freed  from  their  positions  in  the  solid  and  form  a dense  gas. 

The  gas  density  is  the  same  as  that  of  the  solid  for  an  appreciable 
time.  For  breakdown  to  occur,  we  require  that  a chain  reaction  can  take 
place  inside  the  solid  so  that  a gaseous  cavity  is  created  with  the 
dissociation  products.  The  reaction  sequence  as  it  might  occur  in^ 
crystalline  sodium  chloride  will  now  be  described;  the  general  ideas 
are  believed  to  be  applicable  for  all  types  of  chemical  bonding. 
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Within  solid  sodium  chloride,  the  sodium  and  chlorius  atoms  exist 
*4”  •• 

as  ions,  Na  and  Cl  , respectively.  These  might  be  converted  to  gas 

atoms  having  different  states  of  ionization,  such  as  Na°  , Na+  , 

0 + Sas  gas’ 

Cl  , Cl  , and  Cl  , through  the  reactions: 
gas’  gas’  gas’  ° 


1) 

2) 

3) 

4) 

5) 

6) 


Na+  . . . + e -*  Na°  . . , + Q, 
solxd  solid  X1 

Na°  , . . + Q0  -»  Na° 
solxd  2 gas 

+ Q-  -»  Na**"  + e 

gas  3 gas 


Cl  . . . -»  Cl  + Q. 
solxd  gas  4 


Cl  + Qc  -»  Cl  + e 
gas  ^5  gas 

Cl°  + Q.  -»  Cl+  + e. 
gas  ^6  gas 


Reactions  3),  5),  and  6)  occur  in  the  gaseous  state  and  are  activated 
by  collisions  among  the  particles  of  the  gas.  The  entire  reaction 
scheme  must  be  highly  exothermic  if  the  chain  reaction  is  to  be 
vigorous. 


The  above  reaction  scheme  is  depicted  in  a two-dimensional  format 
ii  Figure  1.  Figure  la  shows  a region  with  a single  extra  electron. 

Inis  might  stay  trapped  for  an  appreciable  time,  but  would  not  cause 
lattice  disruption.  In  Figure  lb,  the  electron  density  is  increased 
so  that  the  ionic  charges  on  several  nearest-neighbor  sodium  atoms  are 
neutralized  as  indicated  in  reaction  step  1)  (note  the  five  central 
sodium  atoms  in  the  figure).  Because  of  their  thermal  energies,  the 
neutralized  sodium  atoms  [reaction  step  2)]  can  move  with  ease  from 
their  former  equilibrium  positions  (Figure  1c) . The  neutral  sodium 
atoms  experience  only  slight  binding  energy  to  their  former  equilibrium 
positions.  [Reaction  step  2)  should  occur  more  readily  at  high  tempera- 
tures, suggescing  that  the  breakdown  strength  should  decrease  with 
increasing  temperature.]  The  movement  of  the  sodium  atoms  removes  the 

constraints  on  the  adjacent  Cl  atoms.  Figure  Id  shows  the  system  after 

the  Cl  atoms  have  moved,  through  their  mutual  repulsion,  from  their 
equilibrium  positions  [reaction  step  4)].  In  so  doing  they  lose  some 
potential  energy  and  gain  kinetic  energy;  thus,  by-  subsequent  collisions 
they  heat  up  the  gas  created  by  the  atoms  that  have  left  their  equilib- 
rium posicions.  Figures  le  and  If  show  how  the  electron  density  in  the 
gas  can  be  built  up  through  further  collisions  [reaction  steps  3),  5), 
and  6)].  Finally,  in  Figure  lg,  the  second  stage  of  the  chain  reaction 
is  depicted  after  the  excess  electron  density  has  caused  more  sodium 
atoms  to  be  freed  from  the  solid  matrix. 
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Ig)  REACTION  CONTINUES  WITH 
ADDITIONAL  SODIUM  ATOMS 
ENTERING  GAS. 


Figure  1.  Schematic  description  of  the  bond 
disruption  and  chain  reaction  in  NaCI. 
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Note  that  the  reaction  sequence  has  the  capability  of  releasing 
more  electrons  than  were  needed  for  initiation  and  hence  has  the  poten- 
tial for  a chain  reaction.  The  nature  of  the  mechanism  is  believed  to 
be  quite  general:  nonbonding  should  occur  whenever  the  local  density 

of  electrons  is  made  too  high,  regardless  of  the  detailed  nature  of  the 
chemical,  bond. 


(4)  Establishment  of  a Gaseous  Conducting  Channel.  The 
result  of  the  chain  reaction  is  the  conversion  of  the  solid  in  a local 
region  to  a partially  ionized  gas.  As  the  chain  reaction  proceeds,  a 
gaseous  cavity  is  formed  within  the  solid.  The  reactions  must  take 
place  with  great  rapidity  because  sufficient  energy  must  be  continually 
added  to  the  gas  to  compensate  for  heat  transfer  to  the  surrounding 
medium.  As  the  size  of  the  cavity  increases,  the  medium  about  the 
cavity  is  subjected  to  very  high  stresses  and  cracks  will  start  to  form. 
As  they  form,  new  surface  area  is  exposed  and  the  chain  reaction  spreads 
to  this,  thus  maintaining  the  high  pressure. 

(5)  Completion  of  the  Channel.  We  assume  that  the  break- 
down conduction  current  starts  when  a gaseous  channel  bridges  the  elec- 
trodes. The  currents  prior  to  this  are  associated  with  high  field  pre- 
breakdown conduction  mechanisms;  these  are  assumed  incidental  to  break- 
do™  except  as  they  give  rise  to  a high  local  space  charge  density  that 
can  initiate  the  chain  reaction.  If  the  gaseous  channel  cannot  com- 
pletely bridge  the  electrodes  during  the  time  in  which  the  voltage  is 
applied,  then  the  system  experiences  a "partial  breakdown",  one  in  which 
damage  has  occurred  without  completion  of  the  conducting  path.  If  the 
gaseous  channel  does  bridge  the  electrodes,  conduction  through  the 
channel  begins  since  the  gas  is  ionized.  As  conduction  proceeds,  fur- 
ther erosion  of  the  walls  can  occur  because  of  the  high  electron  density 
and  the  heating  of  the  plasma  by  the  external  circuit.  Experimentally, 
the  conductivity  of  the  system  in  the  breakdown  region  changes  by  a 
factor  of  about  10  [10]  in  a time  of  about  10  nsec,  measured  from  the 
observed  beginning  of  the  gaseous  channels  to  the  onset  of  breakdown 
conduction  [5,30], 

b.  Simplified  Theory  for  the  Growth  of  a Gaseous  Channel 

The  reaction  kinetics  of  real  systems  depend  on  the 
detailed  reactions  present.  As  was  shown  for  the  hypothesized  reactions 
for  breakdown  in  NaCl,  many  reactions  are  possible  during  the  dissocia- 
tion and  recombination  processes  associated  with  creation  of  the  gaseous 
channel.  In  a dielectric  chat  is  chemically  complex,  many  more  dissocia- 
tive and  recombinatorial  processes  would  be  involved.  However,  one  of 
the  significant  features  of  dielectric  breakdown  is  that  the  breakdown 
fields  for  homogeneous  specimens  of  different^  materials  are  remarkably 
similar.  Thus  it  becomes  plausible  to  think  in  terms  cf  a model  having 
highly  simplified  chemical  kinetics.  Such  a model  will  be  outlined  in 
this  section.  The  solutions  for  this  model,  which  depend  on  the  treat- 
ment of  the  heat  flow  probLem  and  detailed  assumptions  on  kinetic 
coefficients  have  not  yet  been  explored. 
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We  assume  that  the  critical  density  of  charge  necessary  for 
initiating  the  chain  reaction  is  present;  the  proposed  theory  is  intended 
to  describe  the  growth  of  the  gaseous  cavity.  The  rate  of  reaction 
wherein  the  solid  is  converted  to  gas  is  assumed  to  depend  upon  the 
electron  density  n£  in  the  gas  and  on  the  wall  area  of  the  cavity.  For 

each  molecule  removed  from  the  solid,  some  neutral  atoms,  positive  ions, 
negative  ions,  and  electrons  are  added  to  the  gas.  Each  material  will 
have  its  own  peculiar  distribution  of  these  reaction  products.  Also, 
there  is  some  total  energy  released,  partially  as  kinetic  energy  of  the 
individual  particles  and  partially  as  electron  excitation,  when  a mole- 
cule of  the  solid  is  converted  to  the  gaseous  phase.  In  our  simplified 
model,  we  describe  the  breakup  without  identifying  the  atomic  species 
present.  When  an  average  molecule  breaks  away  from  the  solid,  we  will 
say  that  a certain  number  of  neutrals,  positive  ions,  negative  ions, 
and  electrons  are  added  to  the  gas.  Atomic  and  ionic  species  will  not 
be  further  differentiated;  i.e.,  an  atom  of  sodium  in  the  gas  will  be 
treated  the  same  as  an  atom  of  chlorine  in  the  gas,  a positive  sodium 
ion  the  same  as  a positive  chlorine  ion,  etc.  For  purposes  of  describing 
the  gas,  we  will  assume  its  mass  density  is  the  same  as  that  of  the  solid 
and  the  average  volume  per  atom  in  the  gas  is  the  same  as  the  average 
volume  (volume  of  a molecule  in  the  solid/No.  of  atoms  in  one  molecule) 
of  an  atom  in  the  solid. 

We  define  wall  generation  coefficients  for  describing  the  dissocia- 
tion at  the  cavity  walls  as  follows: 

1)  = fraction  of  wall  atoms  released  as  positive  ions 

2)  OC  = fraction  of  wall  atoms  released  as  negative  ions 

i 

3)  a = fraction  of  wall  atoms  released  as  neutrals  1 

° 1 

4)  a = fraction  of  the  number  of  wall  atoms  of  electrons  [ 

released  at  the  wall.  \ 

f 

The  wall  generation  coefficients  are  assumed  to  be  independent  of  the  I 

gas  parameters.  For  simplicity  the  ions  are  assumed  to  be  singly  j 

ionized. 

The  gas  (neutrals,  positive  ions,  negative  ions,  and  electrons)  is  [ 

a dynamic  system.  Particles  are  added  continuously  and  the  volume  of  | 

gas  increases  as  the  solid  is  eroded.  At  any  instant  the  cavity  gas  j 

will  be  assumed  to  be  sufficiently  close  to  equilibrium  values  that 
there  is  meaning  to  speaking  of  temperature,  pressure,  and  particle 
concentrations  as  system  parameters.  Since  there  is  heat  flow  from  the 
cavity  gas  while  the  cavity  is  growing,  the  temperature  and  pressure 
will  vary  with  time.  Thus  the  densities  of  neutrals:,  positive  ions, 
negative  ions,  and  electrons  (n.  n+,  n_,  and  n^,  respectively)  will 

vary  through  the  processes  of  recombination,  electron  attachment  and 
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detachment,  and  ionization  within  the  gas.  The  volume  processes  and 
their  kinetic  coefficients  are  as  follows: 

1)  Ionization  of  neutrals,  kinecic  coefficient 

neutral  ->  (+  ion)  + electron 

2)  Recombination,  kinetic  coefficient  7^ 

(+  ion)  + electron  -»  neutral 

3)  Charge  transfer  from  a negative  ion  to  a positive  ion, 
kinetic  coefficient  7 ^ 

(-  ion)  + (+  ion)  2 neutrals 

4)  Charge  exchange  between  two  neutrals,  kinetic  coefficient  7^ 

neutral  + neutral  -»  (+  ion)  + (-  ion) 

5)  Electron  attachment  to  a neutral,  kinetic  coefficient  7^ 

electron  + neutral  -*  (~  ion) 

6)  Electron  detachment  from  a negative  ion,  kinetic 
coefficient 

(-  ion)  ->  electron  + neutral. 

The  kinetic  coefficients  are  functions  of  temperature  and  pressure. 

The  equations  that  determine  the  state  of  the  cavity  system  are 
taken  to  be  the  following: 

1)  Charge  neutrality  condition 

n + n = n , 

2)  Constancy  of  density  of  solid  and  gas 

n . + n + n = a , where  a = number  of  atoms  per 

unit  volume  of  the  solid 

3)  Rate  of  growth  of  the  cavity: 

a)  The  dielectric  is  assumed  isotropic  and  the  cavity 

spherical.  The  electron  density  is  taken  to  be  uniform 
at  all  points  on  the  surface  of  the  cavity.  The  number 
of  particles  added  to  the  cavity  per  unit  time  is 
assumed  proportional  to  the  electron  density  in  the 
cavity  and  the  wall  area;  that  is, 


1- 


d(a  V)/dt  = kn  4rtr2 
„ v ' e 

where  V is  the  volume  of  the  cavity,  r is  the  cavity  radius,  and  lc  is 

a proportionality  constant. 

b)  The  proportionality  constant  is  assumed  to  be  independent 
of  the  gas  parameters.  " n^  is  one  of  the  dependent  vari- 
ables of  the  system  and  its  value  as  a function  of  time 

is  part  of  the  solution.  Thus,  using  V=(4/3)rtr  , the 
above  equation  can  be  integrated  to  obtain  r at  any  time 
t.  The  result  is 

r = (k/av)  J ne(t)  dt 
0 

4)  Rate  equations: 

a)  Positive  ions 

dttj/dt  = a+kne4rtr2  + f^n  + - 7^^  - /2n+n_ 

b)  Negative  ions 

2 2 

dn  /dt  = a_kne4rtr  + 7^nen  + 72n  ~ ^2n-  ” ^2n+n- 

c)  Neutrals 

2 

dn/dt  - CX^kn^Jtr  + 7^r-+ne  + 72n+n  + P2n- 

- w ■ v • v2 

d)  Electrons 

2 

dn  /dt  = a kn  4rtr  + p.n  + (3,n  - y.n.n  - y.nn 
e'  ee  ^2  - K1  '1  + e '4a 

5)  Heat  flow  from  cavity  to  surroundings: 

a)  The  temperature  of  the  cavity  changes  according  to  the 
balance  between  the  energy  added  due  to  the  wall  reac- 
tions and  the  energy  conducting  into  the  surrounding 
medium.  The  heat  capacity  of  the  cavity  gas  becomes 
strongly  temperature-dependent  when  the  temperature  is 
high  enough  to  excite  electrons  and  produce  ionization. 
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kne4rtr“ea  = c(T)  M dT/dt  + dQ/dt  , 

where  is  the  average  energy  released  into  the  gas  per  atom  converted 

to  gas  at  the  walls,  c(l)  is  the  specific  heat  capacity  of  the  gas,  M 
is  the  mass  of  gas  in  the  cavity  at  time  t,  T is  the  temperature  of  the 
cavity  gas,  and  dQ/dt  is  a negative  quantity,  the  heat  flow  from  the 
gas  to  the  surrounding  medium. 


b)  The  heat  flow  from  the  gas  to  the  surroundings  is 

determined  by  solving  the  heat  conduction  equation  for 
the  spherically  symmetric  system  with  the  heat  flow  to 
the  surrounding  medium  equal  to  dQ/dt  in  the  above 
equation.  The  heat  conduction  equation  for  spherical 
symmetry  is 


Si  _ ..  l S 

St'  « ~2  Sr 
r 


c)  The  temperature  at  the  cavity  surface  is  the  same  as  that 
inside  the  cavity.  (A  detailed  solution  using  the  heat 
conduction  equation  is  probably  unattainable,  but  a 
simplified  "law  of  cooling,"  such  as  Newton’s  law  of 
cooling,  might  work  adequately.) 


I I 

II 
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I 
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The  initial  conditions  are: 

at  t = 0,  r = 0,  n+  = 0,  n_  = 0,  n = 0,  ne  = (ne>CRIT  * 

Here,  (ne)GRjq.  is  the  charge  density  necessary  to  initiate  the  chain 
reaction. 

The  method  of  solution  xs  numerical  integration.  First,  reasonable 
values  must  be  chosen  for  the  coefficients.  The  temperature  dependence 
of  the  P’s  and  y’s  and  c(T)  are  required.  The  constraint  equations, 

n + n = n.  and  n,  + n + n = a , 
e - + + v 

can  be  used  to  eliminate  two  of  the  four  dependent  variables  in  two  of 
the  rate  equations.  Thus,  we  might  eliminate  n+  and  n_  in  the  rate 

equations  for  ng  and  n.  The  nQ  and  n rate  equations  would  have  to  be 

solved  simultaneously  with  the  cavity  site  equation  and  the  heat  flow 
equations.  Ultimately,  the  solution  would  give  the  cavity  temperature 
as  a function  of  time,  the  cavity  radius  as  a function  of  time,  n£(t) 

and  n(t).  With  these  quantities  and  using  the  constraint  conditions  in 
the  rate  equations  for  n+  and  n , each  of  these  equations  can  be  expressed 

in  terms  of  one  dependent  variable  and  functions  of  time  and  numerically 
integrated. 
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